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1.  Introduction 


An  ultra-wideband  radar  and  a  plane-wave  antenna  have  been  used  to  measure 
the  high-resolution  scattering  response  of  bare  saline  ice,  snow-covered  ice  and 
pancake  ice  during  the  winter  seasons  of  1994  and  1995  at  the  US  Army  Cold 
Regions  Research  and  Engineering  Laboratory  (CRREL).  The  objectives  of  these 
experiments  were  to  study  various  mechanisms  for  simulating  roughness  and  to 
understand  scattering  mechanisms  better. 

These  backscatter  measurements  were  made  at  2-18  GHz  and  0.5-16.5  GHz 
during  the  1994  and  1995  experiments,  respectively,  and  were  made  for  incidence 
angles  ranging  from  0°  to  55°.  These  broadband  measurements  provide  a  vigorous  test 
for  models  at  many  frequencies  using  a  single  system  instead  of  many  single-frequency 
systems. 

2.  The  Plane-Wave  System 

In  order  to  perform  broadband  scattering  measurements,  we  developed  a 
network-analyzer-based  ultra-wideband  (UWB)  radar  using  a  compact  range  antenna. 
This  system,  as  well  as  providing  exceptional  range  resolution,  also  provides  excellent 
angular  resolution  because  of  the  plane  waves  generated  by  the  compact  range  antenna 
[Gogineni  et  al.,  1995].  Previous  radar  systems  using  hom  antennas  had  beamwidths  in 
the  range  of  6°  to  20°,  which  made  it  difficult  to  determine  the  scattering  mechanisms 
from  distributed  targets  as  a  function  of  incidence  angle.  With  high  range  resolution,  on 
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the  order  of  centimeters,  we  were  able  to  resolve  the  reflections  from  the  complex  ice 
layers,  which  was  never  accomplished  before. 

2.1  Compact  Range  Antenna 

To  obtain  better  resolving  capability  and  angular  response  from  saline  ice  we 
used  a  compact  range  antenna  with  a  wideband  (2-  to  18-GHz)  network-analyzer-based 
radar.  Previous  radar  systems  used  antennas  whose  illumination  of  the  distributed  target 
contained  a  wide  range  of  incidence  angles,  which  made  it  difficult  to  determine 
primary  scattering  sources.  Figure  2.1  shows  a  comparison  between  the  plane-wave 
system  and  conventional  systems. 


(a)  (b) 

Figure  2.1.  (a)  Conventional  radar  systems,  (b)  Plane-wave  system. 


A  42-inch-diameter  parabolic  reflector  antenna  was  used  to  transform  spherical 
waveforms  from  the  offset  AEL  horn  feed  antenna  into  plane  waves  of  uniform  phase 
over  an  area  approximately  equal  to  the  area  of  the  reflector  and  at  ranges  up  to 
0.5*D2/X,  where  D  is  the  reflector  diameter  and  A,  is  the  wavelength.  The  feed  was 
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adjusted  to  place  the  phase  center  of  the  feed  at  the  focal  point  of  the  reflector.  Both 
vertically  and  horizontally  polarized  measurements  can  be  made  with  this  system  by 
properly  orienting  the  feed.  The  rims  of  the  reflector  and  the  feed  support  were  covered 
with  echosorb  to  minimize  the  distortion  of  the  incident  plane  wave  and  to  reduce  the 
scattering  from  these  objects.  Figure  2.2  shows  the  field  of  the  plane  wave.  The  distance 
indicated  in  the  figure  is  relative  to  the  edge  of  the  parabolic  dish. 


I - 1 - 1 - 1 - 1 _ I _ I _ I 

5  10  15  20  25  30  35  40 


Distance  (in.) 

Figure  2.2.  Field  of  plane  wave  6  feet  away  from  the  antenna  shown  for  1995  system. 


2.2  Step-Frequency  Radar 

The  step-frequency  radar  [Izuka  et  al.,  1984]  was  used  together  with  the 
compact  range  antenna  to  measure  the  scattering  coefficient  of  young  saline  ice.  The 
main  feature  of  the  step-frequency  radar  is  the  high  resolution  that  can  be  obtained  by 


3 


exploiting  the  phase  property  of  the  signal.  The  transmitter  signal  is  stepped  in 
frequency  beginning  over  the  desired  bandwidth,  and  the  amplitude  and  phase  of  the 
received  signal  are  measured  at  each  frequency.  The  measured  data  are  concatenated 
and  Fourier-transformed  to  obtain  the  impulse  response  of  the  target.  The  power  return 
as  a  function  of  frequency  is  obtained  by  range  gating  to  isolate  the  return  from  clutter 
and  taking  the  inverse  Fourier  transform. 

2.2.1  Description 


Phase  Lock 


Figure  2.3.  Block  diagram  of  network-analyzer-based  step-frequency  radar. 


The  figure  above  shows  a  simple  block  diagram  of  the  network-analyzer-based 
step-frequency  radar  used  in  our  experiments.  The  synthesized  source  can  generate  a 
swept  signal  in  the  range  of  50  MHz  to  40  GHz.  The  start-and-stop  frequency  can  be 
easily  selected  within  this  range  via  software.  To  achieve  high  frequency  accuracy,  a 
portion  of  the  source  signal  is  routed  to  the  R  sampler  in  the  receiver  and  routed  back  to 
the  source  for  phase  locking.  The  signal  separation  device  consists  of  a  power  splitter, 
two  directional  couplers  and  a  transfer  switch.  The  power  splitter  couples  some  of  the 
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incident  signal  to  the  R  sampler  for  phase  locking  and  reference.  The  transfer  switch  is 
used  to  switch  between  measurements  in  the  forward  direction  (Sn  and  S21)  and  reverse 
direction  (S22  and  S12).  The  directional  couplers  couple  the  reflected  signal  from  the 
target  to  the  A  sampler  in  the  receiver  for  Si  1  measurement  and  to  the  B  sampler  for  S2i 
measurements  when  the  transfer  switch  is  set  in  the  forward  direction.  The  receiver 
then  downconverts  the  received  signal  to  IF. 

In  the  experiments  that  we  conducted,  the  frequency  of  the  network  analyzer 
was  stepped  through  16  GHz  of  bandwidth  (2-18  GHz  in  1994  and  0.5-16.5  GHz  in 
1995)  in  10-MHz  steps  resulting  in  0.94-cm  range  resolution.  Figure  2.4  shows  the 
system  configuration  of  the  ultra-wideband  (UWB)  radar  system.  Scattering 
measurements  were  made  in  Sn  mode  with  this  system  at  an  outdoor  saline  pond  at  the 
CRREL  facility  in  Hanover,  New  Hampshire. 


Figure  2.4.  Setup  of  the  UWB  radar-compact  range  antenna  system. 
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3.  US  Army  Cold  Regions  Research  and  Engineering  Laboratory 


Experiment  (CRRELEX)  Description 

To  understand  the  backscatter  mechanisms  of  sea  ice,  radar  backscatter 
measurements  were  performed  on  artificially  grown  sea  ice  in  an  indoor  facility  at  the 
US  Army  Cold  Regions  Research  and  Engineering  Laboratory  (CRREL)  during  the 
1994  and  1995  winter  seasons.  The  frozen  pond  was  made  with  saline  water  in  order  to 
simulate  sea  ice.  The  pond  is  about  15m  wide  by  20m  long.  The  facility  has  a  catwalk 
on  railroad  tracks  that  were  used  to  place  the  antenna  mount,  the  network  analyzer  and 
the  data  acquisition  system.  The  mobility  of  the  catwalk  permitted  us  to  obtain  four 
independent  samples  at  each  angle  and  also  make  200  independent  measurements  at  a 
fixed  angle  for  probability  density  function  (pdf)  calculations. 

In  January  1994,  measurements  were  made  on  bare  ice  and  snow-covered  ice  at 
incidence  angles  ranging  from  0°  to  55°  with  VV  polarization.  Pdf  measurements  were 
also  made  on  this  ice  at  nadir.  In  March  1994,  measurements  were  made  on  bare  ice 
and  snow-covered  ice  at  incidence  angles  ranging  from  0°  to  60°  with  VV  and  HH 
polarizations.  In  January  1995,  measurements  were  made  on  pancake  ice  at  incidence 
angles  ranging  from  0°  to  50°  with  VV  and  HH  polarizations.  Pdf  measurements  were 
also  made  on  the  pancake  ice  at  0°  and  23°.  To  create  pancake  ice,  waves  were 
generated  on  open  water  with  a  motorized  paddle.  The  wave  action  caused  the 
formation  of  frazil  ice  first,  and  this  slowly  grew  to  pancake  ice  during  freeze-up. 
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The  tables  below  summarize  the  measurements  made  during  these  experiments. 


Table  3.1.  CRRELEX  ’94  January  Data  Summary 


Date 

8  Jan 

0°-55° 

9  Jan 

0°-55° 

9  Jan 

0°-55° 

10  Jan 

0°-55° 

10  Jan 

0°-55° 

11  Jan 

0° 

12  Jan 

0° 

14  Jan 

0° 

14  Jan 

0°-55° 

15  Jan 

0° 

15  Jan 

0°-50° 

15  Jan 

0°-50° 

15  Jan 

0° 

16  Jan 

0° 

E 


DI3  AO  fvj  >10 


133  (Pdf) 


103  (Pdf) 


128  (Pdf) 


95  (Pdf) 


126  (Pdf) 


108  (Pdf) 


ISEjSSHHI 

Polarization 

Snow-covered 

ice 

VV 

2-  18  GHz 

Snow-covered 

ice 

w 

2  -18  GHz 

Bare  ice 

VV 

2 -18  GHz 

Snow-covered 

ice 

VV 

2- 18  GHz 

Bare  ice 

VV 

2- 18  GHz 

Snow-covered 

ice 

w 

2  - 18  GHz 

Bare  ice 

VV 

2  - 18  GHz 

Rough  ice 

w 

2  -18  GHz 

Rough  ice 

VV 

2  -18  GHz 

Snow-covered 

ice 

VV 

2- 18  GHz 

Snow-covered 

ice 

VV 

2- 18  GHz 

Rough  ice 

VV 

2 -18  GHz 

Rough  ice 

VV 

2  - 18  GHz 

Bare  ice 

VV 

2  - 18  GHz 

Table  3.2.  CRRELEX  ’94  March  Data  Summary 


Date 

MUM 

E5S33!29NHi 

Polarization 

11  Mar 

0°  -  60° 

4 

Snow-covered  ice 

VV&HH 

2  -  18  GHz 

12  Mar 

0°  -  60° 

4 

Bare  ice 

VV&HH 

2  - 18  GHz 

Table  3.3.  CRRELEX  ’95  Data  Summary 


Date 

Angles 

M 

Ice  Type 

Pol. 

Frequency 

5  Jan 

0°-50° 

4 

Bare  ice 

VV 

1  - 16  GHz 

7  Jan 

0°-40° 

4 

Snow-covered  ice 

vv 

0.5- 16.5  GHz 

20° 

3 

Snow-covered  ice 

VV 

0.5-  16.5  GHz 

17°,23°, 

31°,41o,44o,50° 

1 

Bare  ice 

vv 

2  - 18  GHz 

10  &  11  Jan 

0°-50° 

3 

Open  water  to 
Pancake  ice 

vv 

0.5- 16.5  GHz 

23  Jan 

0°-50° 

4 

Pancake  ice 

vv 

0.5- 16.5  GHz 

24  Jan 

0° 

200  (Pdf) 

Pancake  ice 

vv 

0.5  - 16.5  GHz 

25  Jan 

0°-50° 

4 

Pancake  ice 

vv 

0.5  - 16.5  GHz 

& 

HH 

25  Jan 

23° 

200  (Pdf) 

Pancake  ice 

VV 

0.5  - 16.5  GHz 

*  Bistatic  Measurements 


8 


4.  Signal  Processing  of  the  Ultra- Wideband  Radar  Data 


The  primary  parameter  to  be  extracted  from  our  measurements  is  the 
backscattering  coefficient  of  the  simulated  sea  ice  as  a  function  of  the  incidence  angle 
and  frequency.  The  radar  equation  for  the  plane-wave  system  is 


qO  _  cal  ^cal 

^ cal  A m 


(4.1) 


where  o°  is  the  backscattering  coefficient  of  the  ice, 

Pr  is  the  power  returned  from  the  ice, 

Pcai  is  the  return  power  from  a  target  of  known  radar  cross  section,  acai, 

Ain  is  the  area  illuminated  by  the  antenna  on  the  target,  and 
Acai  is  the  area  illuminated  by  antenna  on  the  calibration  target. 

This  equation  differs  from  the  conventional  radar  equation  because  of  the 
absence  of  range  dependence  for  the  plane  wave.  The  radar  equation  for  the  plane  wave 
system  is  independent  of  range  because  of  its  narrow  beamwidth.  To  compute  the 
illuminated  area,  we  first  have  to  determine  the  distance  between  the  center  of  the  field 
(similar  to  those  in  figure  2.2)  and  the  point  where  the  magnitude  falls  off  by  6  dB  for 
both  the  horizontal  and  vertical  fields  over  the  frequencies  of  operation.  The  distance 
between  the  6-dB  points  is  the  distance  used  in  the  illuminated  area  calculations  (4.2). 


The  vertical  and  horizontal  fields  have  been  measured  in  the  near-field  region  (6  ft.) 
from  1  to  20  GHz  in  steps  of  1  GHz  for  both  vertical  and  horizontal  polarizations  at  The 
Ohio  State  University  [Nassar,  1992].  The  horizontal  and  the  vertical  distances  in 
meters  are  shown  in  figure  4.1(a-d)  for  the  antennas  used  during  the  1994  and  1995 
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CRREL  experiments.  The  distances  have  been  fitted  with  a  linear  fit  to  interpolate  for 
frequencies  at  which  the  patterns  were  not  measured.  The  vertical  diameter  for  the 
illuminated  area  is  also  modified  by  the  incidence  angle.  Hence,  the  illuminated  area  is 


4*  /(/)  = 


4cos(0) 


(4.2) 


where  Dy(f)  is  the  vertical  diameter, 

Dn(f)  is  the  horizontal  diameter,  and 
0  is  the  incidence  angle. 


Figure  4.1(a).  Horizontal  diameter  used  for  illuminated  area  calculation  in  1994. 
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Figure  4.1(b).  Vertical  diameter  used  for  illuminated  area  calculation  in  1994. 


Figure  4.1(c).  Horizontal  diameter  used  for  illuminated  area  calculation  in  1995. 
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Figure  4.1(d).  Vertical  diameter  used  for  illuminated  area  calculation  in  1995. 


To  extract  the  scattering  coefficient  of  the  ice,  we  need  first  to  determine  the 
return  power  as  a  function  of  frequency.  The  chart  in  figure  4.2  below  illustrates  the 
steps  involved  in  extracting  the  frequency  response  of  the  return  from  ice,  Pr(f). 


Figure  4.2.  Data  processing  steps  to  extract  frequency  response  of  return  from  ice. 


12 


To  determine  the  location  of  our  target,  we  need  first  to  obtain  the  high- 
resolution  spectrum  of  our  data.  In  the  estimation  of  a  spectrum,  there  are  a  variety  of 
windows  that  we  can  choose  to  minimize  leakage.  Unfortunately  most  conventional 
windows  that  have  good  leakage  resistance  usually  suffer  from  poor  resolution.  David 
J.  Thomson  [1982]  introduced  a  technique  of  estimating  high-resolution  spectra  using 
multiple  windows.  These  windows  are  optimum  in  the  sense  of  maximizing  the 
concentration  of  signal  energy  within  a  specified  bandwidth,  W.  The  windows  are 
discrete  prolate  spheroidal  sequences  (DPSS,  also  known  as  Slepian  sequences)  that 
satisfy  a  Toeplitz  matrix  eigenvalue  equation  given  by 


N^1  sin2jtW(n  —  m)  (k) 
m= o  7c(n  -  m ) 


0 N,W)  =  \(N,W)v(nk\N,W ) 


where  N  is  the  number  of  data  points, 

W  is  the  bandwidth,  which  is  in  the  order  of  1/N, 
k  is  the  column  number, 

v®  are  the  eigenvectors  that  correspond  to  the  windows  to  be  used  to  maximize 
the  energy  concentration  within  the  specified  bandwidth,  W,  and 
are  the  eigenvalues  that  give  us  an  estimate  of  the  leakage  that  the  data  set 
incurs  for  the  window,  v®  . 

We  only  need  the  first  K=2NW  terms  with  the  largest  eigenvalues,  but  Thomson 
suggests  that  we  use  the  first  2NW-1  terms  to  minimize  leakage  due  to  the  higher  order 
window  leakage. 
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The  complex  spectrum,  S(t)  (t  is  the  travel  time  to  the  target),  can  be  estimated 
as  follows 

Xk(t)  =  IFFT{x(f)xvw}  (4.4) 

S(0  =  7:Z,7-XJt(0  .  (4.5) 

A  k=  04 

where  x(f)  is  the  raw  step-frequency  data,  and 

Xk(t)  is  the  individual  spectrum  estimate  corresponding  to  the  data  window  v®. 

4.1  Coherent  Noise  Reduction  (CNR) 

Since  the  data  were  collected  in  a  single  antenna  mode,  the  spectrum  will 
include  the  reflection  from  the  antenna  and  feedthrough  (figure  4.3).  These  signals  are 
usually  much  stronger  than  those  of  the  target,  especially  at  higher  incidence  angles 
(>10°).  At  these  angles,  the  return  from  the  ice  is  barely  visible. 

IFFT  of  windowed  and  zero  padded  data 


Range  (m) 

Figure  4.3.  Impulse  response  of  ice. 
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Usually,  there  are  two  ways  that  these  unwanted  reflections  are  reduced.  The 
first  way  is  to  subtract  coherently  the  background  return  from  that  of  the  ice.  This  can 
be  accomplished  by  pointing  the  antenna  to  the  sky,  but  we  could  not  do  this  because  of 
the  limited  extension  range  of  the  actuator.  The  second  way  is  to  implement  a  hardware 
range  gate  to  suppress  these  unwanted  reflections,  but  this  technique  requires  very  fast 
switches  for  short  ranges.  We  tried  range  gating  with  a  commercial  hardware  range 
gate,  the  LINTEK  Range  Gate,  but  we  were  unsuccessful  using  this  hardware.  To 
reduce  the  feedthrough  signals  we  employed  the  coherent  noise  reduction  (CNR) 
technique  [Beaven,  1995],  The  primary  assumption  of  this  technique  is  that  the  phase 
of  the  return  from  a  distributed  target  varies  randomly  from  one  sample  to  another  while 
the  phase  of  the  feedthrough  signals  remains  constant.  CNR  is  carried  out  by  coherently 
averaging  several  independent  samples  at  each  angle  and  subtracting  the  average  from 
each  sample.  Coherent  averaging  reduces  the  return  from  the  target  because  of  the 
varying  phase  and  gives  us  an  estimate  of  the  system  clutter,  which  is  coherent.  If  we 
have  N  independent  samples  (Vi...VN),  CNR  applied  to  Vi  is  simply 

Vicnr  =  -  2^/  (4.6) 

Figure  4.4  shows  an  example  of  CNR  done  on  a  15°  return  from  pancake  ice. 
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CRRELEX  95  -  Return  from  Pancake  Ice  (VV  Pol.) 


Relative  Range  from  Radar  (m) 


Figure  4.4.  Result  of  CNR  on  pancake  ice  at  15°  incidence  angle.  These  data  were 
collected  using  the  plane- wave  system  during  the  CRRELEX  ’95 
experiment. 

4.2  Filtering  &  FFT 

To  isolate  the  ice  return  from  other  signals  in  the  range  domain,  we  needed  to 
design  a  narrow-band  bandpass  digital  filter.  We  designed  the  narrow-band  filter  using 
Thomson’s  Multiple  Windows  to  obtain  the  maximum  concentration  of  energy  within  a 
specified  bandwidth,  W.  The  filter,  H,  is  obtained  by  taking  the  FFT  of  the  windows  for 
the  designed  bandwidth,  W,  dividing  them  by  the  leakage  factor  and  summing  them. 
The  steps  are  as  follows 

Fk  =FFT{vw}  (4.7) 
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(4.8) 


1  1 

H=-iTFk 

K  k=  0  Ak 


The  following  is  the  filtering  sequence  used  to  extract  the  wideband  response  of 

the  return. 

1)  From  the  spectrum  of  the  signal  determine  the  return  associated  with  the  target 
and  synthesize  a  Thomson’s  bandpass  filter  centered  on  the  return  from  the  ice. 

2)  The  FFT  of  the  filtered  signal  will  give  us  the  frequency  response  of  the  ice 
return.  The  drawback  to  this  response  is  the  effect  of  the  filter’s  transient 
response  to  the  range-domain  band  edge,  where  the  transfer-function  response  is 
assumed  to  go  abruptly  to  zero. 

3)  To  reduce  the  ringing  effect,  we  created  a  correction  factor  for  the  frequency 
response  [Higgins,  1990].  This  is  done  by  simulating  the  step-frequency  data 
with  a  reflection  coefficient  of  one  with  the  range  being  set  equal  to  that  of  the 
target.  The  same  process  of  windowing,  zero  padding  and  DFFT  performed  on 
the  measured  data  is  repeated  for  the  simulated  data.  The  filter  used  for  the 
simulated  data  should  be  identical  to  that  of  the  target.  The  FFT  of  this  filtered 
signal  is  the  correction  factor  that  is  divided  by  the  frequency  response  of  the 
target.  Figure  4.5  (a)  and  (b)  shows  the  effects  of  filtering  simulated  data  with 
and  without  the  correction  factor.  The  filter’s  passband  passed  only  4  range 
cells  (4*AR).  The  corrected  frequency  response  has  about  2%  error  at  the  edges. 
Notice  the  ringing  effect  in  the  beginning  of  the  uncorrected  filtered  response. 
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-‘Actual  response;  --'Uncorrected filtered  response 


Figure  4.5(a).  Frequency  response  of  filtered  data  using  Thomson’s  filter.  ( —  =  actual 
response;  —  =  uncorrected  filtered  response). 


-'Actual  response;  --'Corrected  filtered  response 


Figure  4.5(b).  Filtered  data  with  correction  factor  applied  ( —  =  actual  response; 
—  =  corrected  filtered  response). 
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The  effectiveness  of  the  correction  is  dependent  on  the  width  of  the  filter 
passband.  The  larger  the  filter  passband,  the  more  accurate  will  be  our  frequency 
response,  but  widening  the  filter  may  include  spurious  responses  from  our  measured 
data.  For  low  incidence  angles  (0°-10°),  two  to  four  range  cells  in  the  passband  are 
sufficient  to  obtain  the  frequency  response.  For  higher  incidence  angles,  a  wider 
passband  should  be  used  to  include  the  returns  from  the  larger  illuminated  area.  The 
corrected  frequency  response  of  the  target,  Pr(f),  is  used  in  (4.1)  to  compute  the 
backscattering  coefficient.  Figure  4.6  shows  an  example  of  the  frequency  response, 
Pr(f),  of  the  return  from  pancake  ice  at  0°. 


Figure  4.6.  Uncalibrated  return  power  as  a  function  of  frequency  for  pancake  ice  at  0°. 
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5.  CRRELEX  Results 


The  results  of  the  CRREL  experiments  are  presented  in  three  parts.  The  first 
part  consists  of  a  few  impulse  response  plots  of  the  various  ice  types;  namely,  bare  ice 
and  snow-covered  ice  at  a  few  angles.  The  second  part  consists  of  angular  responses  of 
the  scattering  coefficients  for  bare  ice,  snow-covered  ice  and  pancake  ice  types  at  C 
band  and  Ku  band.  A  comparison  between  the  scattering  coefficients  obtained  using  the 
plane-wave  and  the  spherical-wave  antennas  are  shown  here.  Comparison  with  field 
measurements  for  pancake  ice  is  also  made  in  this  section.  The  third  part  consists  of 
frequency  responses  of  the  scattering  coefficients  at  various  incidence  angles.  A 
comparison  between  the  measured  and  theoretical  scattering  response  of  bare  ice  at  5°, 
10°  ,15°  ,20°  and  30°  is  also  shown  here. 

5.1  Impulse  Response 

Figure  5.1  shows  the  typical  response  from  bare  ice  at  nadir.  The  return  from 
smooth  bare  ice  at  nadir  resembles  a  delta  function.  There  is  a  spurious  response  8  cm 
after  the  ice  reflection.  This  could  be  due  to  the  system  response  since  a  similar 
response  is  also  observed  from  the  snow-covered  ice.  As  the  incidence  angle  is 
increased  to  5°,  the  return  power  drops  by  almost  30  dB  (figure  5.2).  This  demonstrates 
the  excellent  angular  resolution  that  we  are  able  to  achieve  with  our  system. 
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CRREL  94  -  Return  from  Bare  Ice  (W  pol.) 


Figure  5.1.  Return  from  bare  ice  at  0°. 


CRREL  94  -  Return  from  Bare  Ice  (W  pol.)  @  0  &  5  degrees 


Figure  5.2.  Return  from  bare  ice  at  0°  and  5°. 


Figure  5.3  demonstrates  the  excellent  resolving  capability  of  the  plane- wave 
system.  Notice  the  clear  reflections  from  the  air/snow  and  the  snow/ice  interfaces, 
which  are  about  5  cm  apart.  The  reflection  from  the  snow/ice  interface  is  about  seven 
times  higher  than  that  from  the  air/snow  interface.  This  is  due  to  the  higher  dielectric 
contrast  between  snow  and  saline  ice  than  between  snow  and  air. 
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CRREL  94  -  Return  from  Snow-Covered  Ice  (W  pol.) 


Figure  5.3.  Return  from  snow-covered  ice  at  0°. 


We  used  a  resonant  probe  [Smith  and  King,  1974]  to  measure  the  dielectric 
constants  of  ice  and  snow.  The  resonant  frequency  of  the  probe  was  measured  in  free 
space  and  the  medium.  The  relationship  between  the  resonant  frequency  of  the  probe  in 
the  medium  and  its  dielectric  constant  is  given  by 

^med 

£  . 

c'a\r 

where  emcd  &  £air  are  the  dielectric  constants  of  the  medium  and  air,  respectively,  and 
fmed  &  f  air  are  the  resonant  frequencies  of  the  probe  in  the  medium  and  air,  respectively. 

The  dielectric  constants  of  ice  and  snow  were  measured  to  be  3.5  and  1.17, 
respectively.  The  relationship  between  the  reflection  coefficient  and  the  dielectric 
constant  of  a  material  is  given  by 


f  •  ' 

J  air 

i  f  ,  i 


(5.1) 
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(5.2) 


r_  yi75~-i 

where  T  is  the  reflection  coefficient  at  the  dielectric  interface  and  £1  &  £2  are  the 
dielectric  constant  of  medium  1  and  medium  2,  respectively. 

The  reflection  coefficients  at  the  snow/ice  interface  and  air/snow  interface  were 
found  to  be  0.266  and  0.0406,  respectively.  The  ratio  of  the  reflection  coefficient  from 
the  ice  to  snow  is  6.56:1.  This  agrees  with  our  measured  result  shown  in  figure  5.3. 

5.2  Angular  Scattering  Response 

The  scattering  responses  were  computed  and  calibrated  to  a  metal  sphere  with  a 
radar  cross  section,  acai,  of  -13.9  dBm2  over  the  entire  frequency  range. 

Figure  5.4  illustrates  how  wide-beam  (spherical-wave)  antennas  corrupt  the 
scattering  response  from  bare  ice  at  near  vertical  incidence.  The  scattering  coefficient 
obtained  using  the  plane-wave  system  at  5°  is  25  dB  lower  than  that  obtained  using  the 
conventional  system.  The  results  of  the  hom  antenna  are  corrupted  by  scattering  from 
vertical  until  the  incidence  angle  increases  to  about  20°.  The  data  collected  using  the 
hom  antenna  needs  to  be  corrected  for  the  large  beamwidth  in  order  to  obtain  the  correct 
angular  response.  This  result  confirms  the  superiority  of  the  plane-wave  system  over 
conventional  systems  for  measuring  the  angular  dependence  of  the  scattering  response 
for  vertical  incidence. 
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Figure  5.4.  Comparison  of  backscattering  coefficient  from  bare  ice  obtained  using 
plane-wave  and  spherical-wave  antennas  at  Ku  band. 

Figures  5.5  and  5.6  below  show  the  angular  response  of  several  ice  types  at  C 
band  and  Ku  band,  respectively.  From  the  plots  above  we  can  see  that  the  scattering  for 
pancake  ice  is  higher  than  bare  ice  for  angles  greater  than  zero.  The  results  show  that 
there  is  more  incoherent  scattering  from  pancake  ice  because  of  its  rougher  surface.  The 
scattering  from  snow-covered  ice  is  slightly  higher  than  that  from  bare  ice  because  of 
the  slightly  rougher  surface.  We  can  also  see  for  the  three  ice  types  that  surface 
scattering  may  dominate,  at  least  until  30°,  because  of  the  scattering  coefficient  falloff 
with  angle. 
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Figure  5.5.  Angular  response  of  bare  ice,  snow-covered  ice  and  pancake  ice  at  C  band 
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Figure  5.6.  Angular  response  of  bare  ice,  snow-covered  ice  and  pancake  ice  at  Ku 
band. 
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Figure  5.7  shows  the  agreement  between  field  [Beaven,  1992]  and  laboratory 
backscatter  characteristics  of  pancake  ice  at  35°. 


Incidence  Angle  (degrees) 


Figure  5.7.  Pancake  ice  comparison  at  C  Band. 


5.3  Frequency  Scattering  Response 

At  0°  (figure  5.8),  we  see  that  scattering  from  pancake  ice  is  increasing  with 
frequency.  The  scattering  level  increased  from  21-29  dB  as  the  frequency  was  increased 
from  2-15  GHz.  When  the  frequency  is  increased,  the  surface  seen  by  the  radar 
becomes  rougher.  This  would  mean  that  the  scattering  in  the  specular  direction  will 
begin  to  reduce  and  the  scattering  will  begin  to  increase  in  other  directions.  The 
scattering  in  other  directions  are  called  the  incoherent  components.  From  figure  5.8,  we 
can  infer  that  at  the  initial  frequency  of  2  GHz,  the  incoherent  components  have  been 
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dominating  for  the  pancake  ice  and,  as  the  frequency  is  increased,  we  get  an  increase  in 


scattering  from  the  incoherent  components.  As  for  the  snow-covered  ice,  it  can  be 
observed  that  the  scattering  decreased  from  25-19  dB  as  the  frequency  was  increased 
from  3.5-16.5  GHz.  From  this  result  we  can  infer  that  the  coherent  component  has  been 
dominating  at  the  initial  frequency  of  3.5  GHz  and  as  the  frequency  is  increased,  we  get 
a  decrease  in  scattering  from  the  coherent  component  and  an  increase  in  scattering  from 
the  incoherent  components.  The  higher  scattering  at  the  initial  frequencies  could  be  due 
to  wet  snow  at  the  snow-ice  interface.  As  for  the  bare  ice,  there  is  no  significant  change 
in  the  scattering  level  across  the  frequency  range,  which  means  that  there  is  only  a 
dominant  coherent  component.  To  explain  the  relative  scattering  levels  of  these  three 
ice  types,  more  modeling  effort  is  required  to  determine  the  contribution  from  the 
incoherent  and  coherent  components. 
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Figure  5.8.  Frequency  response  of  bare  ice,  snow-covered  ice  and  pancake  ice  at  0°. 

When  the  incidence  angle  is  increased  to  5°  (figure  5.9),  there  is  a  15-  to  20-dB 
drop  in  the  scattering  coefficient  of  pancake  ice  and  about  a  25-  to  30-dB  drop  in  the 
scattering  coefficient  of  bare  ice  and  snow-covered  ice.  The  large  drop  in  the  scattering 
level  is  due  to  the  reduction  in  coherent  scattering.  The  high  scattering  coefficient  at  the 
initial  frequency  is  due  to  the  coherent  contribution  of  the  ice.  This  is  due  to  the  larger 
beamwidth  of  the  radiating  pattern  at  the  lower  frequencies.  The  scattering  level  for 
pancake  ice  increased  as  a  function  of  frequency.  This  is  an  indication  of  dominant 
incoherent  scattering. 
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Frequency  response  of  Bare,  Snow  cover  and  Pancake  ice  @  5° 


Figure  5.9.  Frequency  response  of  bare  ice,  snow-covered  ice  and  pancake  ice  at  5°. 

The  results  for  the  bare  ice  measurements  at  5°  and  beyond  were  compared  to 
theory.  The  coherent  component  was  computed  using  Fung  and  Eom’s  [1983]  model 
and  the  incoherent  component  was  computed  using  the  integral  equation  method  (IEM) 
model  [Fung,  1994].  The  solution  to  determine  the  coherent  contribution,  a°c ,  is  given 
by  Fung  and  Eom  as 

r 

°c  =  ^rexp(-4fc  2  a2 ) exp(-02  /  B 2 )  (5.3) 

where  Tp  is  the  reflection  coefficient  for  polarization  p  (either  H  or  V), 

0  is  the  incidence  angle, 
k  is  the  wavenumber, 
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Gh  is  the  rms  roughness,  and 


B  is  given  by 


B2  = 


(  R\2 


(khfiy 


where  (3  is  the  one-sided  beamwidth  and 
h  is  the  range  to  the  scatterer. 

The  incoherent  component,  G°nc ,  using  the  IEM  model  is  given  as 


o'.  =  — ex 

me  2 


where  kz  =  kcosO,  kx=ksin0, 


C  =  (2kzoffvv  exp ,(-fcz2cr2)  +  M  M  M  + 


(5.4) 


(5.5) 


(5.6) 


and  the  symbol  W(n)(-2kx,0)  is  the  Fourier  transform  of  the  nth  power  of  the  surface 
correlation  function.  Other  symbols  are  given  in  Appendix  2A  of  Fung,  1994. 

The  theoretical  scattering  coefficient,  o°v,  is  given  by 


_  O  o  .  _  o 

°vv  —  ®inc 


(5.7) 


The  parameters  used  for  these  models  are  given  in  Table  5. 1 . 
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Table  5.1.  Model  parameters  for  bare  saline  ice 


Parameter 

Value 

rms  roughness,  Gh 

1.0  mm 

correlation  length,  Lc 

2.2  cm 

dielectric  constant,  £r 

3.38-j0.22 

One-sided  beamwidth  @  3 
GHz 

2.3° 

Equivalent  spherical  source 
range,  h 

9.3  m 

The  beamwidth  and  range  are  factors  only  for  the  coherent  contribution  case. 
The  equivalent  spherical  source  range  was  calculated  using  the  radius  of  the  illuminated 
area  from  figure  4.1(a-d)  and  the  beamwidth.  From  figure  5.10,  we  can  determine  h 
using  the  trignometric  relation 

k  ~  Tan(p)  (5'8' 


< - > 

2r 

Figure  5.10.  Range  calculation  for  equivalent  spherical  source. 
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To  determine  the  beamwidth  of  the  radiating  pattern,  we  measured  the  impulse 
response  of  the  plane-wave  system  using  a  mesh  screen  at  0°,  0.5°,  1.0°,  1.5°,  2.0°  and 
5.0°.  The  power  return  from  the  mesh  screen  at  angles  beyond  0°  is  due  to  the  coherent 
component.  From  the  coherent  returns  at  the  various  incidence  angles,  we  can  estimate 
the  radiation  pattern  of  the  antenna.  Figure  5.11  shows  the  one-sided,  two-way  radiating 
pattern  of  the  plane-wave  antenna  at  3  GHz,  which  has  been  fitted  with  a  Gaussian 
curve.  The  one-sided  beamwidth  for  the  one-way  pattern  is  the  -6dB  point  of  the  two- 
way  radiating  pattern.  At  3  GHz,  this  corresponds  to  about  2.3°. 


Figure  5.11.  Radiating  pattern  of  the  plane-wave  antenna  at  3  GHz. 


Figure  5.12  shows  the  comparison  between  the  theoretical  scattering  values  and 
the  measured  values  of  bare  saline  ice  at  5°.  From  this  figure  we  see  that  the  coherent 
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component  is  dominating  the  scattering  source  until  about  6  GHz  and  the  incoherent 
scattering  begins  to  dominate  beyond  6  GHz. 


Figure  5.12.  Comparison  between  theoretical  and  measured  scattering  of  bare  ice  at  5°. 

At  30°  (figure  5.13),  the  scattering  level  for  bare  ice  and  snow-covered  ice 
increased  by  about  12  dB  with  increasing  frequency.  This  increase  corresponds  to  an 
increase  in  incoherent  scattering.  The  scattering  from  pancake  ice  also  exhibits  a  similar 
behavior  but  with  a  higher  scattering  level  due  to  its  rougher  surface. 

Figure  5.14  shows  the  scattering  comparison  between  the  theoretical  scattering 
values  and  the  measured  values  of  bare  saline  ice  at  30°. 
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Figure  5.13.  Frequency  response  of  bare  ice,  snow-covered  ice  and  pancake  ice  at  30‘ 
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Figure  5.14.  Comparison  between  theoretical  and  measured  scattering  of  bare  ice  at 


Figures  5.15-5.17  show  the  comparison  between  the  theoretical  and  measured  values  of 
bare  ice  at  10°,  15°  and  20°.  We  believe  that  the  discrepancies  at  the  initial  frequencies 
at  these  angles  are  due  to  the  coherent  component  from  the  sidelobes  of  the  radiating 
pattern. 


Figure  5.15.  Comparison  between  theoretical  and  measured  scattering  of  bare  ice  at 
10°. 
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Figure  5.16.  Comparison  between  theoretical  and  measured  scattering  of  bare  ice  at 
15° 


Figure  5.17.  Comparison  between  theoretical  and  measured  scattering  of  bare  ice  at 
20°. 
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6.  Conclusions  and  Future  Work 


An  ultra-wideband  radar  and  a  plane-wave  antenna  have  been  used  to  measure 
the  high-resolution  scattering  response  of  bare  saline  ice,  snow-covered  ice  and  pancake 
ice  at  the  US  Army  Cold  Regions  Research  and  Engineering  Laboratory  (CRREL) 
during  the  1994  and  1995  winter  seasons. 

We  have  developed  signal  processing  techniques  using  Thomson’s  Multiple 
Windows  to  estimate  the  high-resolution  spectrum  and  to  extract  the  wideband  response 
of  the  signal.  We  have  shown  how  the  ringing  effect  associated  with  narrowband  filters 
can  be  reduced  by  using  a  correction  factor. 

The  results  of  our  measurements  show  the  excellent  spatial  and  angular 
resolution  that  can  be  achieved  with  the  plane-wave  system.  The  scattering  coefficients 
obtained  using  the  plane-wave  system  show  that  surface  scattering  dominates  at  least 
until  30°.  At  0°,  the  scattering  from  bare  ice  stayed  fairly  constant  with  about  1-  to  2-dB 
variation  across  the  2-  to  18-GHz  frequency  range,  while  the  scattering  from  pancake 
ice  increased  from  21  to  29  dB  with  increasing  frequency.  On  the  other  hand,  the 
scattering  from  snow-covered  ice  decreased  from  25  to  19  dB  with  increasing 
frequency.  This  shows  that  there  is  more  contribution  from  incoherent  components  for 
pancake  ice  at  higher  frequencies,  while  for  snow-covered  ice  there  is  less  contribution 
from  the  coherent  component  at  the  higher  frequencies.  At  5°  incidence  angle,  we  can 
see  that  the  coherent  component  is  still  dominating  up  to  6  GHz,  after  which  the 
incoherent  component  begins  to  dominate.  At  30°  incidence  angle,  the  scattering 
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increased  as  a  function  of  frequency,  which  is  an  indication  of  dominant  incoherent 
scattering.  With  increasing  incidence  angle,  the  scattering  level  dropped.  This  is  an 
indication  of  dominant  surface  scattering.  Field  measurement  of  pancake  ice  at  35° 
shows  agreement  with  laboratory  measurements. 

We  need  to  ascertain  the  radiating  pattern  of  the  antenna  beyond  5°  to  determine 
if  there  are  any  coherent  sources  of  scattering  at  these  angles.  We  also  need  to  do 
modeling  on  the  different  ice  types  to  determine  the  relative  surface  and  volume 
contributions.  There  is  also  a  potential  for  using  a  focussed  beam  for  better  penetration 
and  probing  of  surface  roughness  over  the  ice. 
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